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ABSTRACT 

We present the first metallicity distribution functions of the old field populations in the Magellanic 
Clouds. Our metallicities are based on the Fourier decomposition of Typeab RRLyrae light curves 
from the Optical Gravitational Lensing Experiment (OGLE-HI). On the metallicity scale of Zinn & 
West; we find a mean metallicity of [Fe/H] = —1.50 ± 0.24 dex based on 16776 RRLyrae stars in the 
Large Magellanic Cloud (LMC). For the Smah Magellanic Cloud (SMC) we obtain -1.70 ± 0.27 dex 
based on 1831 RRLyrae stars. These uncertainties represent the intrinsic spread in the population 
rather than the standard deviation of the mean. 

Our results are in good agreement with the few existing spectroscopic metallicity determinations 
for LMC RRLyrae stars from the literature. For both the LMC and the SMC the metallicity spread 
exceeds 1 dex in [Fe/H]. The distribution of metallicities in both Clouds is very uniform, and no 
significant metallicity gradient is detectable. We also do not find any pronounced populations of 
extremely metal-poor RRLyrae candidates with metallicities well below —2 dex, although we need to 
caution that the photometric method used may overestimate the metallicities of metal-deficient stars. 
Moreover, because of stellar evolutionary effects one does not expect to observe many RRLyrae stars 
among very metal-poor horizontal branch stars. We suggest that the Magellanic Clouds experienced 
fairly rapid and efficient early enrichment involving pre-enriched gas as well as possibly gas infall, while 
metal loss through outflows does not seem to have played a significant role. Moreover we suggest that 
the differences in the metallicities of the old population of LMC and SMC make an origin from a 
single, common progenitor unlikely, unless the separation happened very early on. 
Subject headings: (Galaxies:) Magellanic Clouds — Stars: abundances - Stars: variables: RRLyrae 



1. INTRODUCTION 

The origin and nature of the Magellanic Clouds (MCs) 
pose an interesting puzzle. They have long been known 
as the closest irregular satellites of the Milky Way, al- 
though galaxy surveys reveal that star-forming satel- 
lites in close proximit y to massive galaxies are rare (e.g., 
while very accurate Hubble Space 



James & Ivory 



2011) 



Piatek et al. 



(e.g. 



the MCs dKallivayalil et al.||2006b[ 
have evoked a new discussion of their orbits 
|et al. 2007, 2010). Wit h recent new measurem ents of 



20081 
., 'Besla 



lower p roper motions by Costa et al. ( |2011| ) and [Vieira 
et alj (2010) the discussion has been further intensi 
In 



tied. In particular, these studies have raised questions of 
whether the MCs mi ght be on their first passage around 
the Milky Way (e.^.,[L uxet al.|20"TQ|and re ferences cited 
above; but see also, e.g., Diaz &; Bekki 2011 ). Similarly, it 
remains unclear whether the MCs were formed together, 
whether the Small Magellanic Cloud (SMC) might have 
been captured by the Large Magellanic Cloud (LMC) 
later on, or whether interactions or mergers might be 
responsible for the ch emistry and peculia r structure of 
the SMC (e.R.JMathe ws Qn et al. 1988 ; Kallivayalil et al. 
2006a[ [Costa et al.||2000HTsujimoto fc Bekki 2009)7 



Une of the possible avenues of exploring the early 
history of the MCs relies on the study of their old popu- 
lations. Generally, both globular clusters as well as field 
stars are suitable tracers of old populations. The LMC 
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does indeed contain a number of old globular clusters 
that are similarly old as the oldest globulars in other 
Galactic satellites and in the Milky Way (e .g., Brocato] 
etal]p96| [Olsen et al.||T998l [Grebel fc Galla gher 200jr 



In terms of their chemical composition, the LMC's gloBu- 
lar clusters span a similar range of metallicities and show 
simila r [g/Fe] ratios as stars in the Galactic halo (e.g.. 



Telescope- bdised determinations of the proper motions of l 



Johnson et al 



2006 



Mucciarelli et al. 



2010) 



trast, the SMC only contains one old globu. 



In con- 
ar cluster, 

NGC121, which is a few Gyr youn ger than the oldes t 
Galactic and LMC globular clusters ( [Glatt et al.|2008b[ ). 
Thus one would ideally wish to turn to more numer- 
ous tracers of old populations such as those potentially 
present among field stars. However, genuinely old field 
stars are difficult to identify in galaxies that experienced 
long- lasting star formation such as the MCs (e.g., Holtz- 



man et al. 



1999 Gahart et al. 20081 Sabbi et al. 



iff 



since stellar populations of many diiierent ages are then 
superimposed, rendering otherwise popular stellar trac- 
ers such as red gi ants rather ambigu ous indicators (see. 



e.g.. Figure 10 in Koch et al. 2006). Using star clus 



ters with spectro scopic metalli c ities an d ages from main 



sequence fitting, Glatt et al. (2008a) showed that the 
SMC exhibits a considerable spread in metallicity at any 
given age. The study of evolutionary histories is further 
exacerbated if the galaxy in question has a non-negligible 
depth extent as appears to be the case for the SMC. Un- 
der such circumstances, using easily identifiable, special 
types of star s as tracers is a useful approach (see, e.g., 
Grebel|1997| Fi gureS). One of these evolutionary tracers 



2 



Haschke et al. 



for evaluating the parameters of the old field population 
are RRLyrae variable s. With ages of at least 9 Gyr (e.g. 
Sarajedini et al.|2QQ6 ), they allow us to explore the early 
stages of the evolutionary history of a galaxy. 

When well-defined light curves are available, RRLyrae 
stars can easily be identified. Apart from their distance 
and spatial distribution, their metallicity is a parameter 
of particularly high interest. Either spectroscopic mea- 
surements or Fourier decomposition of their li ght curves 



can be used to deduce the met allicity [Fe/H] ( jKovacs fc 
' Zsoldos|p 95l lJu rcsik fc Km^cs 1996 ). 



Old populations including KK Lyrae stars tend to be 
widely distributed and often form the most extended 
stellar populations in galaxies. This provides us with 
the possibility to not only measure the mean metallicity 
of genuinely old stars but also to obtain the metallicity 
spread, to derive detailed metallicity distribution func- 
tions (MDFs), and to identify possible metallicity gradi- 
ents. Overall, such studies can provide a fairly detailed 
picture of the chemical evolution that a galaxy experi- 
en ced at early times. 

Butler et ah] ([1982) first used the AS method ( |Pre- 



ston| |l959|) to investigate the metallicity of six RRLyrae 



stars m the LMC and ten RRLyrae stars in the SMC. 
They found a significant difference between the two sam- 
ples, with a mean metallicity of [Fe/H] = — 1.4 ± 0.1 dex 
for the LMC and -1.8 ± 0.2 dex for the SMC. In more 
recent studies, spectroscopic metallicities of old stars in 
the LMC were de termined for sets of 74 to 98 RRLyrae 



stars (Gratton et al.||2004l |Borissova et al.||2004l 12006 



see TablejT]), while lor the SMC no further spectroscopic 
data are available. These spectroscopic studies, which 
are concentrated close to the bar of the LMC, revealed 
mean metallicities of about [Fe/H] = —1.50 dex. 

With microlensing surveys such as the Optical Grav 



itational Lensing Experiment (OGLE) ( [Udalski et al. 



1992[ |1997{ |2008D, carefully chosen, extended regions 
5tella 



with high stellar densities have been regularly monitored. 
In the Magellanic Clouds, OGLE focussed on the dense 
central regions. These surveys are well suited for finding 
variable stars, such as RRLyrae stars. Thanks to the 
large number of repeat observations, variability proper- 
ties of these stars can be determined with high precision. 
The parameters of the observed light curves of RR Lyrae 
stars of type ab ca n be used to determine th e metallicity 
as show n first by Kovacs fc Zsoldos (1995). [Jurcsik fc' 
Kovacs (|l996) presented a formalism to derive metallici- 
ties as a Function of the period and the Fourier phase ^31. 
For cases where the Fourier parameters c anno t be deter- 



mined , for in stance Brown et al. (2004) and Sarajedini 
et al. (2006) presented equations to infer the metallic- 



ity using only the period. These approaches have larger 
uncerta inties than the equations presented by [Jurcsik fc] 



Kovacs 



( 1996} and are not used in our study. 
In recent years, the public data of the OGLE survey 
have been used to determine the MDF and the metal- 
licity gradient of the MCs with RRLyrae stars. The 
OGLE II data were used by 'Smolec' f2005) to investi- 
gate the metallicity distribution of 5451 LMC RRLyrae 
stars and the differences in [Fe/H] for Blazhko and non- 
Blazhko stars. For the non-Blazhko RRLyrae stars a 
mean metallicity of [Fe/H] = —1.218 ± 0.004 dex was de- 
termined, while the Blazhko stars had a mean value of 
[Fe/H] = -1.28 ±0.01 dex. ,Deb fc Singh (2010) used the 



OGLE II data of the SMC to determine metallicities for 
335 RR Lyrae stars of type ab and found a mean metal- 



licity of [Fe/H] = -1.56 ± 0.25 dex. In Feast et al. (2010 



from now on FAWIO) a gradient of 0.01 ±0.002 dex kpc 
was inferred for the RRLyrae stars of the LMC, observed 
by OGLE HI. 

We use the entire sample of RR Lyrae stars provided 
by the third data release of the OGLE collaboration to 
calculate individual metallicities for each star of the MCs. 
We discuss the origin of our data in Section|2j In Sec- 
tion|3]we present the method for calculating metallicities 
usedin this paper. The results for the LMC are given in 
Section!!] and for the SMC in Section[5l We thus obtain 
an MDF of the old population and investigate the exis- 
tence/absence of a metallicity gradient for both MCs. We 
test whether the photometric metallicities lead to similar 
results as spectroscopic values in Section[6] and compare 
different methods to estimate the metallicity from pho- 
tometric data in Section[7l Section[8] discusses the results 
obtained. 

2. DATA 

In 1992 the OGLE experiment started its first phase of 
observations of the MCs using the 1 m Swope telescope o f 
Las Campanas Observatory, Chile ( [Udal ski et al.| |1992 ). 
This telescope was replaced by a 1.3 m telescope owned 
by the OGLE collaboration at Las Campanas in 1995, 
with which the second phase of the OGLE experiment 
was performed. The CCD camera of the first two phases 
had a field of view of 15' x 15' and 2048 x 2048 pixels 



(Udalski et al.||1997|). Data were taken preferentially in 
the I band and additionally in the V and B bands un- 
til 2001, when the third phase of OGLE was started. 
The eight CCDs of the OGLE HI camera, each with 
2048 X 4096 pixels, took data until 2009. With a to- 
tal of 8192 X 8192 pixels, t his camera covered a total area 
of 35' X 35' (Udalski 2003). In the LMC nearly 40 square 
degrees were monitored, while in the SMC 14 square de- 
grees were observed using V and / band filters. Alto- 
gether photometric data of high precision are provided 
for 35 million stars in the LMC and 6.2 million stars in 
the SMCR In addition to the full photometric catalogue, 
subsets of interesting objects such as RRLyrae stars or 
Cepheids have also been compiled by the OGLE collab- 
oration. 

In the LMC the OGLE collaboration identified 17693 



RRLyrae stars of typeab (Soszyhski et al. 2008), while 
in the SM C only 1933 RRL yrae stars of typeab were 
found (Soszyhski et al. '2010'). As well as light curve pa- 
rameters, mean magnitudes in the V and / bands and 
periods were determined by OGLE. The / band light 
curves were decomposed using a fifth order Fourier se- 
ries following Simon & Clement (1993) and the Fourier 
amplitudes and phases determined] 

3. METHOD OF METALLICITY DETERMINATIONS 

Determining metallicities of a sample of this size is cur- 
rently only feasible via photometric estimates. We there- 
fore use the techniques of metallicity estimates based on 
Fourier-decomposed light curve parameters. 



Simon & 



^ The catalogues are available from |http : //ogle . astrouw. edu. 
Pl/| 
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Lee ( |1981 ) introduced this method for variable stars, us- 
ing the tirst few orders of a Fourier decomposition, 



N 

V = Ao ^^Ak sm{kujt - 

/c=l 



(1) 



where V is the observed magnitude, Aq^j. is the ampli- 
tude, the phase, uj the angular frequency, and t the 
specific time at which a data point was measured. N is 
the order to which the Fourier series is developed. 

The Fourier parameters of the different orders of 
the decomposition contain information about important 
physical properties of the variables. While the ampli- 
tude Ak is a measure of the skewness, the phase c/)/^ cor- 
responds t o the acuteness of t he light curve ( Stellingwerf| 
fc Dickens 1987). Simon & Clement ( 1993j showed how 
the i:^Qurier'°parameters =0 3 — 3^1 are related to tem- 
perature and luminosity, while Kovacs & Zsoldos (1995) 
presented the close connection between the period, the 
Fourier parameter and the metallicity of RRLyrae 
stars. The relation for a sinusoidal Fourier decomposi- 
tion of a y band l ight c urve is given by Equation (3) of 
Jurcsik fc Kovacs] ( |1996] ). 

l^burier parameters differ depending on the filter in 
which the light curve is observed. As the OGLE pho- 
tometry is mostly done in the / ba nd, the V band 
equatio n by Jurcsik & Kovacs (1996) cannot be used. 
Smolec ( 2005) recalibrated the relation to be suitable for 
I band photometry. Combining light curve parameters 
of 28 RRLyrae stars with their spectroscopic metallici- 
ties ranging from — 1.7 1dex < | Fe/H] TQr; < +0.08 dex, on 
the m etallicity scale of Jurcsik (T995, from now on J95), 
[Smolec (2005) found the linear equation 

[Fe/H]j,, = (-3.142 ±0.646)- (2) 
(4.902 ± 0.375) P + (0.824 ± 0.104) ^31, 

with 031 being the Fourier phase and P being the period 
of the RRLyrae star. We adopt a slightly larger metallic- 
ity range of — 2.0dex < [Fe/H] < +0.2 dex as a presum- 
ably reliable domain. This larger range is the metallicity 
range given by the calibration in Smolec (2005) plus the 
mean uncertainty of the method (±0.18dex). 

The p95 . metalli city scale, which is the basis for the 
relation by | Smolec] (^05), le ads to higher metalli cities 
compare d wit h other scales. [Gratton et al. (2004) find 



that the |J95 | scale is more metal- rich by U.28dex than 
the scale of Z inn fc West' (1984', lat er on ZW84), in 
agreem ent with SanHage (2004) and ^ Clement ini et al. 



( 2005 ) . Using six globular clusters wit h metallicities be- 



tween — 2.25dex < [Fe/H] < — l.Odex, Papadakis et al 



( 2000 ) derived the transformation relation 



[Fe/H] 



ZW84 



1.028 [Fe/H]j93 - 0.242 



(3) 



between the two metallicity scales. Several other metal- 
licity scal es are also co mmonly used. A widely-used scale 
is that of [Carretta fc Gratton (1997). With a new and 
larger sample, Carretta et al. (2009) recahbrated the 
metallicity scale of Carretta fc Gratton (1997) and pro- 
vide a quadratic relation betw een thei r hew metallicity 
scale and the metallicity scale of ZW84 Another popular 



metallicity sc ale was introdu c ed by Harris (1996), which, 
according to Gratton et al. (20041), is indistinguishable 
from .ZW84l 



-1.5 -1 




[Fe/H] 



J95 

Figure 1. Metallicity distribution of LMC RRLyrae stars, com- 
puted based on the Fourier decomposition of their light curves. 
We find a mean metallicity of [Fe/H]^^^^ = — 1.22dex on the scale 
of J95 The bottom axis shows the scale of J95 while the top 
axis shows the scale by ZW84 for which we find [Fe/H] = 
— 1.49dex, using the transformation of Papad akis et al.] ( |2000| ). 
The distribution has a standard deviation ot U.26dex. We also tind 
a few candidates for very metal-poor stars. The error bar shows 
the typical uncertainty of the metallicities for individual stars. 

4. METALLICITIES OF THE LMC 

For all RR Lyrae stars of type ab in the LMC for which 
031 values are available from the OGLE III catalogue, we 
compute metallicities using Equation (pL This results in 
a dataset of 16949 stellar metallicities (Figurejl]). 



4.1. 



Mean metallicities and metallicity distribution 
function 



J95 



scale, we find a mean metallicity 
1.22 dex with a standard deviation of 



Using the 
of [Fe/H] = 

0.26 dex. This translates to a mean metallicity of 
[Fe/H] = -1.50 dex on the ZW84 metallicity scale. The 
uncertainties of the Fourier parameters are not men- 
tioned by the OGLE collaboration. Taking the uncer- 
tainty of the observed magnitude as the 1 a standard 
deviation of its assumed Gaussian error distribution, 
we randomly vary the observed magnitudes to perform 
Monte Carlo simulations of the light curve for each indi- 
vidual star and redetermine the Fourier parameters. We 
find very good agreement between the parameter values 
obtained by OGLE and via our Monte Carlo approach. 
The resulting range of parameters allows us to calculate 
the uncertainty for ^31 and thus of the corresponding 
metallicity of each star. For the mean uncertainty of 
the metallicity we find crintrinsic ~ 0.10 dex. This un- 
certainty is quadratically added to the sy stematic uncer- 
tainty of cTgys = 0.18 dex mentioned by | Smolec (2005) 



and results in a total uncertainty of o-^f^/h] = U.2l dex. 

We convolve the individual metallicity values obtained 
with the derived constant (Gaussian) uncertainty of 
0.21 dex and sum up the resulting Gaussian distributions 
to obtain the "observed" MDF. The solid blue line in Fig- 
ure|2] shows the MDF normalized to 1 (i.e., a probability 
distribution function) . 

Motivated by the near- Gaussian profile of the ob- 
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['=«/"]zW84 
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Table 1 
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Figure 2. Probability distribution function for the metallicity of 
LMC RRLyrae stars. The sohd blue curve shows the sum of indi- 
vidual stellar metallicities convolved with Gaussian uncertainties, 
while the red dash-dotted curve is the likelihood distribution. The 
computation of the metallicities is based on the Fourier decom- 
position of their light curves. We find a mean metallicity of the 
likelihood distribution of [Fe/H]^^^^^ = — 1.23dex on the scale of 
|J95| (bottom axis). The top axis shows the ZW84 scale, with a 
mean value of [Fe/H]^^^^ = — l.SOdex, using the transformation 
of Pa padakis et alT] ( |2000[ ). The MDF has an intrinsic spread of 
0.24 del^^ 

served MDF, we also use a maximum-likelihood approach 
to constrain the underlying MDF of the LMC RRLyrae 
population. Using such an approach allows us to deter- 
mine the intrinsic width of the MDF (assuming a Gaus- 
sian profile) while simultaneously accounting for the con- 
tribution of the measurement uncertainties to the width 
of the observed distribution. The likelihood of our ob- 
served data set originating from an underlying Gaussian 
MDF with its mean, /i, and standard deviation, cr, is 
given by the function 



>c(Ai,cT)=n: 



N 



= 1 727r(a2+a,2) 



exp I 



(-^fef) (4) 



where N is the total number of stars, Xi is the metal- 
licity of star i and cr^ is the associated measurement un- 
certainty. We seek to find the parameters /i and a that 
maximize this likelihood function. 

The dash-dotted line in Figure[2] represents the best 
fit of /i and a and reveals the underlying true MDF with 
a mean metallicity of [Fe/H] = — 1.23dex and an intrin- 
sic width of the distribution of 0.24. 

The highest metallicity found in our sample is 
Fe/Hl^^^ = 2.0 dex, while the lowest value is 
Fe/H]^.^ = —3.81 dex. As explained in Sectional we 
only consider the reliable domain of the calibration and 
thus exclude 107 stars (0.6% of the complete sample) 
at the low-metallicity end and 66 stars (0.4% of the 
complete sample) at the high-metallicity end from the 
dataset when applying this selection. This leaves us with 
a sample of 16776 RRLyrae stars. However, the mean 
and the standard deviation of the metallicity of the re- 
duced sample remain unchanged in comparison with the 
full sample. 

Results on different metallicity scales cannot easily 
be compared. We transform all metallicities to the scale 
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Gratton et al.|(|2004| 
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Figure 3. Distribution of metallicity in the LMC. Metal-poor 
regions are coded in blue, metal-rich regions are red. Overall the 
contour plot is smooth and no gradient or distinct feature is visible. 
For the contour plot, only stars with —2.3 < [Fe/H]2^34 < 0.0 are 
taken into account. The black lines in the middle represent the bar. 
W hile the white asterisk mar ks the optical center of the LMC found 
by |de Vaucouleurs &: Freeman| ( |1972|) , t he black aster isk show s the 
center ot the KHLyrae stars as tound bylHaschke et al. ( submitted! 
|to AJt . 



of[ZW84 using Equation (^3l. This transformation leads 
to a mean metallicity of ^e/H]^gg^^y2W84 ~ —l.SOdex 
on the ZW84 scale, which is in excellent agreement with 
mean values of spectrosco pic measurements of RR Lyrae 
metallicitie s in th e LMC ( Gratton et al.|[2"QQ4 Borissova 
eFaT] [20041 [20061 who used the |Z\V^ nieHllicity scale 



as wellJT 

4.2. Spatial variations in metallicity 

Figure|3] shows the spatial distribution of the metal- 
licities. Only s tars wi th a metallicity from the reliable 
domain of the Smolec[ |2005 ) calibration are used in this 
figure. Individual stars are used for the contour plot, to 
which a grid with constant boxwidth of 1.5' x 1.5' was 
then applied. Figure[3] shows the median metallicity val- 
ues of these boxes. 

In general, a smooth distribution with small-scale vari- 
ations, but without any obvious gradients or distinct fea- 
tures, is visible. Towards the center of the LMC, the den- 
sity of RRLyrae stars increases. Differences in metallic- 
ity are therefore observed on smaller scales in the central 
parts than in the outskirts. This leads to the higher res- 
olution in the center of the figure. The stars that are ex- 
cluded due to their metallicity being above or below the 
adopted validity range are distributed randomly across 
the whole field of OGLE III. 

Figure 4] examines the mean metallicity as a func- 
tion of the distance from the center of the RR Lyrae stars 
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Figure 4. The contours show the total number of RRLyrae stars 
of a given metalhcity as a function of the projected distance rela- 
tive to the center of the LMC as found from RRLyrae stars. The 
mean metallicity does not change, within the errors, with increas- 
ing distance from the center. The grey dots represent the median 
metallicity in bins of 0.5 kpc. 

of the LMC at a = 5^26^ a nd 6 = -69° 75' as found 
bylHaschke et al.| dsubmitted to AJ|). We investigate the 
metahicity and distance distributions in bins of 0.1 dex 
and rings of varying radius. The radii are chosen such 
that all annuli have the same area. The density of stars 
in each box is determined. The contour plot is smoothed 
with a Gaussian filter using 3x3 boxes and a width of 
one box. Additionally we calculate the median metallic- 
ity of each distance bin of 0.5 kpc width, this "running 
median" is represented by the grey dots in Figure|4j 

To determine the gradient of the metallicity, the me- 
dian metallicity is computed in bins of 0.1 kpc and fitted 
using a first-order polynomial. The resulting gradient 
for the whole distribution of —0.03 ± 0.07 dex kpc~^ is 
consistent with zero. If we take only the innermost 8 kpc 
from the center of the RRLyrae stars into account, a gra- 
dient of —0.010 ± 0.014 dex kpc~^ is found. The median 
metallicity seems to be constant in the inner parts of the 
LMC and dropping slightly towards the outskirts. 

To test for non-radial gradients, we divide the 
OGLE III field into 15 similarly-sized fields. As shown 
in Figure |5l the mean metallicity of these fields drops in 
accordance with the very small metallicity gradient ob- 
tained above. The outer fields always show slightly lower 
metallicities than the innermost fields. The median val- 
ue s are nea rly identical to the mean values. 

|FAW10| determined the metallicity gradient of the old 
population using OGLE III RRLyrae stars. These au- 
thors calculated photometric metallicities with two sets 
of equations, both with the period as the only free pa- 
rameter. The uncertainty of this method is therefore 
higher (a = 0.46 dex) than with the method of Jurcsik 



& Kovacs (1996) (a = 0.16 dex). We rec alculate our 



metallicities usmg the equations utilized by [FAWIO and 
compare t hem wit h our estimates (see Table[2| . Equa- 
tion^) of FAWIO and our estimates transformed to the 
|ZW84 metallicity scale lead to similar results when av- 
eraging over all metallicities. For subsets of the sam- 
ple, however, we find considerable differences. We divide 
the sample into subsets where all stars are more metal 
rich than our average or more metal poor. The more 
metal-rich stars have in general higher metallicities in 
our Fourier-based estimates. The metallicities below the 
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Figure 5. The variance of metallicity between the inner- and out- 
ermost fields is tested by dividing the LMC field into 15 similarly 
sized fields. The mean metallicity values on the ZW84-scale are 
shown in black, grey values represent the scale of Jurcsik ( 1995|. 
We see a small decline of the metallicities in the outermost regions 
sampled by OGLE III. The boundaries of the bar are illustrated 
with the dashed grey lines, while the black asterisk marks the cen- 
ter as found by RRLyrae stars and the grey asterisk the optical 
center of the LMC. 

mean are in general lower in our calculations than the 
values computed with Equation (2) by FAWIO, When 
only considering stars at the extremes of ouFlnetallic- 
ity estimates this trend becomes even stronger (Table[2|. 
However, the spatial distribution of the RRLyrae stars 
in each metallicity bin is very similar to the overall dis- 
tribution. 

FAWIO conclude that there is a very small gradient 
of -0.0104 ± 0.0021 dex kpc-^ for the RRLyrae stars of 
the innermost 6 kpc of the LMC. This result is in very 
good agreement with the gradient found in our study. 
Our error estimate is a bit more conservative, as we take 
the intrinsic uncertainties of the metallicity into account. 
Within the errors, we are not able to distinguish between 
no gradient and a very small gradient. 

The result of essentially no metall i city g r adient con- 
firms t he findings of Grocholski et al. (2006); Cole et al. 
(2009) and Sharma et al. ( 2010 )7^ho used either star 
clusters with ages between 1 and 10 Gyr or red giants to 
determine the metallicity gradient of the LMC. A small 
gradient of —0.047 ± 0.003 dex kpc~^ was obtained by 
^iona(2009) using AGB stars, which span an age range of 
3 — 10 Gyr. Based on spectroscopy of red giants 'Carrera 



et al. ([2011) find that the metallicity of these stars in the 
mner okpc of the LMC is roughly constant, but then de- 
creases with increasing radii. This result is supported by 
our findings with the (on-average) older RRLyrae stars. 

5. METALLICITIES OF THE SMC 

In the SMC, Fourier parameters were determined for 
1864 R RLyrae stars of typea b by the OGLE III collabo- 
ration ( Soszyhski et al.||2010[ ). The resulting metallicity 
distribution is shown in l^lgure[6] 

Following the same maximum-likelihood approach as 
described for the LMC RRLyrae, we find the underly- 
ing metallicity distribution function of SMC RRLyrae 
to have a mean value of [Fe/H] = —1.42 dex on the scale 
of J95 with a standard deviation of 0.27 dex as shown 
in Figure[7l As explained in Section|3j we choose the 
metallicity regime between —2.0 < [Fe/H] < +0.2 on 
the scale of IJ95I as the reliable domain of the calibra- 



6 



Haschke et al. 



Table 2 

Differences between metallicity estimates by |FAW10| an d by us f or different ranges of metallicity. The first column gives the metahicity 
scale we adopt, the second column the equation used by |FAW10] in order to derive the metallicity, the third column the metallicity range 
used for the comparison is shown, and the fourth column shows the difference found by calculating A^f = F^/H]hgdjii ~ F^/^FAWio- 



[Fe/H]HOD.Tii scale 


[Fe/H] 


FAWio equation 


[Fe/H] range 


Ahf 


J95 


[Fe/H] = 


-5.62 log 




- 2.81 


-2.00 < [Fe/H] < +0.20 


0.17 


J95 


[Fe/H] = 


-7.82 log 


;P 


-3.43 


-2.00 < [Fe/H] < +0.20 


0.25 


ZW84 


[Fe/H] = 


-5.62 log 


;P 


- 2.81 


-2.00 < [Fe/H] < +0.20 


-0.11 


ZW84 


[Fe/H] = 


-7.82 log 


^.P 


-3.43 


-2.00 < [Fe/H] < +0.20 


-0.03 


ZW84 


[Fe/H] = 


-7.82 log 


\P 


-3.43 


-1.23 < [Fe/H] < +0.20 


0.23 


ZW84 


[Fe/H] = 


-7.82 log 


;P 


-3.43 


-2.00 < [Fe/H] < -1.23 


-0.18 


ZW84 


[Fe/H] = 


-7.82 log 


;P 


-3.43 


-0.50 < [Fe/H] < +0.00 


0.52 


ZW84 


[Fe/H] = 


-7.82 log 


^.P 


-3.43 


-1.00 < [Fe/H] < -0.50 


0.41 


ZW84 


[Fe/H] = 


-7.82 log 


^.P 


-3.43 


-1.50 < [Fe/H] < -1.00 


-0.05 


ZW84 


[Fe/H] = 


-7.82 log 


;P 


-3.43 


-2.00 < [Fe/H] < -1.50 


-0.17 



[Fe/H] 



ZW84 




[Fe/H] 



J95 

Figure 6. The metallicity distribution for the old population of 
the SMC. For 1864 RRLyrae stars [Fe/H] is determined and a 
mean metallicity of [Fe/H]^^^^ = — 1.42dex on the scale of J95 
is obtained (bottom axis), while on the scale of ZW84 we find 
[Fe/H]^g^^ = — 1.70dex (top axis). The standard deviation of 
the distribution is 0.33 dex. The error bar represents the typical 
uncertainty of the metallicities for individual stars. 



tion by lSmolec (2005). All values outside of this calibra- 
tion range are rejected. 27 stars with metallicities below 
[Fe/H] = —2.0 dex are found, while the lowest value is 
measured to be [Fe/H]^. = —3.3 dex. Only 6 RRLyrae 
with metallicities above [Fe/H] = +0.2 dex are present 
and a maximum of [Fe/H]^^^ = 1.7 dex is found. Elim- 
inating these stars, which are randomly distributed in 
space, leaves us with a sample of 1831 RRLyrae stars 
with presumably reliable metallicities. 

The resulting mean metallicity of the old population 
of the SMC is [Fe/H]^_^ = -1.42 dex on the scale of 



J95]and [Fe/H]^^^^ = -1.70 dex on the ZW84 metal- 
licity scale. So far no spectroscopic measurements of 
RRLyrae field stars have been published. We therefore 
compare our result with the metallicity of NGC 121, the 
only globular cluster of the S MC in the a^e ran^e brack- 
eted by our RRLyrae stars. Glatt et al.' ('2008b") found 
an age of 10.5 Gyr for NGC 121 using main-sequence fit 
ting with Dartmou th (Dotter et al. 2007) isochrones. Da 
[Costa fc Hatzidim itriou (1998) used the CaH triplet to 
determine the mean metallicity of NG C 121 an d found 
[Fe/H] = -1.46 ± 0.10 on the scale of |ZW84[ in very 
good agreement with the results we obtain for the old 



[Fe/H] 



ZW84 




[Fe/H] 



J95 

Figure 7. Metallicity probability distribution functions for the 
SMC. The solid blue curve represents the measured metallicities 
convolved with Gaussian uncertainties. A likelihood distribution is 
determined from the solid curve and given by the red dash-dotted 
line. We use 1831 RRLyrae stars and find amean metallicity 
of [Fe/H]^g^^ = —1.42 dex on the scale by |J95| (bottom axis). 
[Fe/H]^g^^ = -1.70 dex is obtained for the scale of |ZW84^ (top 
axis). The intrinsic spread of the distribution is 0.27 dex. 

population o f the SMC 



Recently Kapakos et al. (|2011|) used a sample of 



100 RRLyrae stars from the OGLE III catalogue to de 
termine ph otometric metallicities using the relation by 
Jurcsik & Kovacs (1996). They decomposed the V band 
light curves to determine the Fourier parameters of each 
star. On the metallicity scale of J95 they found a mean 
metallicity of [Fe/H]^^^^ = -1.51 ± 0.41 dex. This is 
again in very good agreement with the metallicities found 
with our sample. 

We use our metallicity estimates on the scale of |J95| 
for the individual RRLyrae stars to plot a contour map 
of the metallicity distribution of the old population (Fig- 
ure|8|. The overall distribution of the metallicity is quite 
smooth and no distinct features are visible. The black 
lines indicate the borders of the OGLE III field. In the 
white regions within this polygon no RR Lyrae stars were 
detected and therefore no contours are drawn. 

The metallicity gradient is measured with respect 
to the center of the SMC determined by K-, M- and 
faint carbon stars (a = 0^51^ and S = -73°^, Gq^ 



nidakis et al. 2009[ ) This center is indicated in FigurejS 
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-1\ 



[Fe/H] 



ZW84 

Figure 8. Distribution of metallicity in the SMC. Metal-poor 
regions are coded in blue, metal-rich regions are red. Overall the 
contour plot is smooth and no gradient or distinct feature is visible. 
For the contour plot only stars with —2.0 < [Fe/H] < +0.2 are 
taken into account. The w hite asterisk represents the cente r found 
by Gonida kis et al.| ( |2009| >, the grey asterisk the center of [Piatek] 
[et al.| ( ,2UUg| r 




Distance [kpc] 



Figure 9. The total number of stars in boxes of metallicity versus 
the projected distance from the center as determined by Gonidakis 
|et al. (2009) (white star in Figure|8j. The grey dots represent the 
median metallicity of small distance bins of 0.5 kpc. Within the 
uncertainties no radial gradient of the metallicity is found. 

with a white asterisk. The spatial density distributi on of 
RRLyrae stars is bimodal ( Haschke et al.|pii prep. ) and 
the determination of a single center is therefore ra ther 
difficult. We have c hosen to adopt the values by |Go- 
[nidakis et al.| (|2QQ9[), because they represent the mean 
center ot the KR Lyr, 



LLyrae distribution quite well. The stel- 
lar density is evaluated in boxes of 0.1 dex with a chang- 
ing radius to keep the total area of each annulus constant. 
The resulting contours are smoothed with a Gaussian 
kernel and shown in Figure|9] 

In order to investigate the presence of a gradient we 
fit a ffist-order polynomial to the running median of the 
SMC. We find a gradient of 0.00 ± 0.06 dex kpc-^ We 
subdivide the OGLE III field of the SMC into 9 fields of 
similar area and determine the mean metallicity in each 
(Figure 10). We note that the mean is very similar to the 
median value. The metallicity stays basically constant, 




ZW84 
J95 



a [hours] 



Figure 10. The SMC field is divided into nine similarly sized 
fields and the mean metallicity for each field is obtained. The mean 
values on the J95 and ZW84 scales are shown in grey and black, 
respectively. The values of the different fields are very similar to 
each other. The black and grey ast erisk indicates the center of the 
SMC as found by |Gonidakis et al.| ( [2009t and Piatek et aL] ( |2008t , 
respectively. 

as expected froni the essentially zero gradient. 

Piatek et al. (12008) used HST proper motion data 
to determine the kinematical center of the SMC. They 
found a = 0^52^8^ and 6 = -72°30^ which is shown 
by the grey star in Figure 8] We do not find pronounced 
radial trends in the metallicity when using this center of 
the SMC either. 

Metallicity gradients of the S MC have been studied 
using s tar clu sters and red giants. |Da Costa fc: Hatzidim- 



itriou 



( 1998 ) used spectroscopic measurements of giants 

and star clusters to search for a poss ible meta llicity gra- 
dient and did not find any evidence. Glatt et al. (2008a) 
studied six intermediate age clusters and found a spread 
of 0.6 dex between the most metal-poor and most metal- 
rich clusters. The spectroscopic metallicity estimates 
lead to the conclusion that no smooth, monotonic age 
metallicity relation is present, but that inst ead there is 
a spr ead in me tallicity at any given age (Kayser et al. 



2007[ Glatt et al. 2008a [). Us ing spectroscopic data o: 



red giants Farisi et al.| ([2009) did not find a gradient 



while Carrera et al. "i 20b8[ ) suggest that one exists using 
a different sampl e oi spectroscopic obs ervations of red 
giant field stars. |Carrera et al.| (|2008) argue that this 
gradient is the result of different ages of the red giants at 
different locations of the SMC. They conclude that fur- 
ther star formation has taken place in the inner parts and 
led to the metallicity gradient among these intermediate- 
age stars with ages between 3 and 10 Gyr. However this 
would not have affected our old RRLyrae stars. 



6. COMPARISON OF SPECTROSCOPIC AND 
PHOTOMETRIC METALLICITIES 

For some RRLyrae stars located in the bar of the 
LMC, spectroscopic metallicities were determined in ear- 
lier studies. 149 LM C RRLyrae wit h spectrosc opic 
metal l icities from Gratton et al .l ([2004) and Boris sova| 
et al.| ( |2004| |2006| ) also have UGLE light curves and are 



1 1 shows the photometric 



thus part oi our sample. Figure 
metallicities we determine with the Fourier decomposi- 
tion method versus the difference of photometric minus 
spectroscopic metallicities from the literature. 
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Figure 11. Comparison of our photometric metallicities with 
spectroscopic metalhcities from the hterature for LMC RRLyrae 
stars. Ah metalhcities are on the |ZW84| scale. Calculating 
A[Fe/H] = [Fe/H]^;^^^ - [Fe/H]^pg^ we find a median difference 
of 0.05 dex for stars with photometric metallicities and their coun- 
terparts in the samples of Gratton et al.| ( ,2004| ) and 
( 2004 2006). The trends m each sample, except for 
(2UU4^, are fitted with a linear polynomial. 



Borissova et al. 



BonssovaeFaT 



Th e spectroscopic measurement s of [Borissova et al.] 
(2004) are based on the method of Layden (1994), who 
used the comparison of the equivalent width ot (JallK 
line with hydrogen Unes. The resu lting metalli city is on 



the scale of |ZW84[ [Borissova eFal | ( |2006j ) and [Gratton 
et al.[ (|2004| used the method ot line mdices (jFreston 



T959) as outlined in Gratton et al. (2004). This method 



Is linked to the metallicity scale oi Harris (1996), which 
has a very small offset of + 0.02 ± 0.01 dex compared to 
ZW84| ( [Gratton et al.|2 QQ?). Due to the nearly negligible 
ottset we do not adjust the metallicity values when ap- 
plying the comparison to the photometric metallicities. 

In FigurepT] we see that there appears to be a sys- 
tematic trend in the deviation between spectroscopic and 
photometric measurements suc h that f or low metallicities 
(lower than —1.8 dex on the | ZW84| scale) the spectro- 
scopic metallicities are typically higher, while for photo- 
metric metallicities higher than ~ —1.3 dex the spectro- 
scopic values are on average increasingly lower. These 



trends are fitted with linear polynomials (Figure 11 and 
Table[3|. The median deviation, when subtractmg the 
spectroscopic metallicity from the photometric metallic- 
ity, is +0.05 dex, while we find a standard deviation of 



0.35 dex. However, the spectroscopic values in the three 
studies cited above show very different characteristics, as 
presented in Table[3l For the 58 RRLyrae from the sam- 
ple by Gratton et al. (2004) and the 28 stars in common 
with Boris sova et al.| ([2004 ) we find a median deviation 
of 0.01 dex and —0.05 dex with a standard deviation of 
0.25 dex and 0.49 dex, respectively, whereas we find a me- 
dian difference of 0.15 dex and a standard deviation of 



0.35 de x for the 66 RRLyrae taken from Borissova et al. 
(2006). They quote an uncertainty range of 0.20 dex for 
their spectroscopic metallicities. 



In [Kunder fc Chaboyer] (|2008 ) the same spectroscopic 
metallicities as used here are compared with the pho- 
tometric data of the V band lightcurves of the LMC 
RRLyrae stars from the MACHO survey. They do not 
find any evidence for a systematic discrepancy between 
the spectroscopic values and the photometric metallic- 
ities c alculated with the relation by Jurcsik fc Kovacsj 
( 1996 our Equation[5|. This leads to the conclusion that 
the systematic diffe rence in our st udy is introduced by us- 
ing the relation by Smolec (2005). However, most of the 
stars are located in the photometric metallicity regime 
where the differences to the spectroscopic measurements 
are quite small. 

We also searched the spe ctroscopically studied, pre- 
sumably old, s tar clu sters of Colucci et ah] ( 2011); Gro-| 



cholski et al. 



^ ([2006") and Muc ciarelh et al.| (|^UUJj) tor 

RRLyrae stars detected by OGLE. We do not find any 
coincidences of RR Lyrae stars in our sample within the 
clusters' diameter. 



7. METALLICITIES FROM TRANSFORMED (/)3i VALUES 



Jurcsik fc Kovacs (1996) used V band light curves 
to determine a relation between (\)'^\ and metallic ity. 



Deb fc Singh 



(12010] 
tH^( 



es- 
:orm 



Dorfi fc Feuchtinger (1999) and^ 

tablished two independent linear relations o 

= + /^^fci transform the / band Fourier param- 
eters to V band values. The values for a and /3 shown in 
Table[4lare taken from Table 3 of Deb fc Singh (2010|. 

^De'B'fc Singh (2010) used the observed light curves of 
the RRLyrae stars of two globular cluste rs, with metal- 
liciti es of -2.16 dex (M68,lLee et al.'2005') and -1.39 dex 
(M3, [Cohen fc Meiendez 2005), to calculate Fourier pa- 
rameters. With these data they obtained values for ol and 
^ that are significantly different from )Dorfi fc Feuchtingerj 
(1999), who used Fourier-decomposed theoretical light 
curve models of RRLyrae stars to determine the rela- 
tion. These models were matched with observational 
data from three globular clusters with metallicities of 
-2.16±0.02dex (M68, Lee et al. 2005), -1.52 ±0.12 dex 
(IC4499, Hankey fc Cole 2011) and -1.18 ± 0.02 dex 
(NGC1851, Carretta et al.,2010) . 

We transform the l^burier parameters provided by 
OGLE using the relations listed in Tableffl The recalcu- 
lated values of (/^)^are inserted into the formula deter- 



mined by [Jurcsik fc Kovacs (1996): 



[Fe/H] = -5.038 - 5.394P + 1.345(^; 



(5) 



to redetermine the MDF of the MCs^ All metallicities 

presented in this section are on the |J95[ scale. Owing 
to the very small number of calibration clusters and the 
narrow metallicity range of Dorfi fc Feuchtingerj ( 1999 ) 
and , Deb fc Singh^(2010J , we use this relation with a note 
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Table 3 

Differences between the photometric and spectroscopic metahicities for a sample of 149 RRLyrae stars. We calculate the mean, median 
and standard deviation for the whole set of stars and for every study. To test for systematics we fit the distribution with a linear 

least-square fit. 





all stars 


Gratton et al. 


(2004) 


Borissova et al. ( 2004 ) 


Borissova et al. ( 2006]) 


mean 


0.037 




-0.008 


-0.041 


0.108 


median 


0.051 




0.011 


-0.052 


0.145 


std 


0.346 




0.254 


0.492 


0.347 


fit 


0.8160[Fe/H] + 1.2427 


0.6181[Fe/H] + 


0.9388 




1.0272[Fe/H] + 1.5884 



Table 4 

Parameters to transform the Fourier coefficient (f)si from the / 
band to the V band. 



= + /3(/>^i 


a 


(3 


|Dorfi^&: Feuchtinger" ('1999') 
peb &; Singh (2010) 


-0.039 
0.436 


0.788 
0.568 



of caution outside the metalhcity regime defined by the 
globular clusters. 

7.1. LMC 



The relation by Dorfi & Feuchtinger| (1999) leads to 



a very different rnetallicity distribution (upper panel 



in Figure 12) than obtained in FigurefT] The resulting 
metallicity distribution is much broader with a stan- 
dard deviation of 0.38 dex, while the mean metallicity 
of [Fe/H]^g^^^j3pgg = —1.20 ± 0.38 dex is very similar to 

the value obtained with t he [Smolecl (|2005 |) rela tion. 

Using coefficients from Deb ^ Singh (2010) to calcu- 
late metahicities of all RR Lyrae stars results in a metal- 
licity distribution (lower panel in Figure 12) similar in 
shape to that shown in Figure[l] The mean metallicity is 
comparable within the errors. With [Fe/Hj^^g^^yj^g-^Q = 

— 1.33 ± 0.29 dex, we find a metal l icity that is about 
0.1 dex lower than from the Smolecj ( |2005| method. For 
this distribution a standard deviation of 0.29 dex is ob- 
tained assuming a Gaussian distribution. 



7.2. SMC 



Using the Dorfi fc Feuchtinger (1999) relation, the 
mean metallicity, on the J95 scale, is [l^e/Hj^g^^^j^pgg = 
— 1.43 dex with a standard deviation of 0.40 dex. As with 
the LMC, the shape of the SMC metallicity distribu- 
tion is much fi atter and wider than with the relation by 
Smolec| ( |2005 >), while the mean metallicity is very similar 
(upper panel of Figureplsl)^ 



The relation of Deb'~E' Singh (2010) reveals a similar 
shape of the metallicity d istribution com pared with that 
based on the relation by Smolec (20051). However, the 
metallicity is about 0.1 dex lower with [I'e/Hlmean/DSio ~ 
— 1.53 dex, while the standard deviation is 0.32 dex. 
These results are shown in the lower panel of Figure[T3] 
and are analogous to the outcome for the LMC. 

Overall, 72 % of the SMC metalli cities determined with 
the relation by Deb & Singh' ('2010') and 21% of the metal- 
licities from the relation by Dorfi &; Feuchtinger (1999) lie 
outside of the metallicity range defined by the calibrat- 
ing chisters used to derive the relations. The relation 
by Smolec (2005|) has not been calibrated for such low 
metallicities as the ones we find for some of our RR Lyrae 
stars, but includes a wider metallicity range. Especially 



[Fe/H],w84 

-1.5 =1_ 




[Fe/H] 
[Fe/H] 



'J95 
'zW84 




[Fe/H] 



'J95 



Figure 12. The upper panel shows the histogram of the LMC 
me talhcity distribution o btained using the conversion of cf)^^ to (f)Yi 
by |Dorfi Feuchtinger] jl999 J, w hile the lower panel is obtained 
witn tne relation ot |ueb & aingn| ( [20To] > . 

in the metallicity interval of —1.71 dex < [Fe/Hjjg^ < 
— 1.30 dex, whe re most of the RR Lyrae stars are located 
in our sample, |Smolec (2005) has a good cover age with 
13 stars, thus justifying our use of Smolec (2005) relation 
as opposed to the other relations described here. 



8. CONCLUSIONS 

We derive photometric metallicities for 16949 
RRLyrae stars in the observed OGLE III field of the 
LMC using the Fourier decomposition of the light curves 
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Figure 13. iDorfi &: Feuchtinger| |T999t conversion of the Fourier 
parameter ^g^^ to (f)^-^^ is used in the upper panel to determine the 
metalhcity of the SMC RRLyrae stars. In the lower panel the rela- 
tion by Deb &; Singh (2010) is applied to determine the metallicity 
distribution. 



( Jurcsik & Kovacs'1996). On the J95| metallicity scale we 
lind a mean metallicity of [Fe/H]^ = — 1.23dex. This 
corres ponds t o a metallicity of [Fe/Hj^^g^^ = — l.SOdex 
on the |ZW8"4| scale, which is in excellent agreement with 
the spectroscopic measurements of LMC RRLyrae stars 



-jv rae stars 

in the literature (Borissova et al. 2004 |2006{ iGrattonJ 
|et al. 2004). Our results are a lso in excellent agreement 
with those of |Smolec| ( |2005| ). The mean metallicity 
of his dataset^ which is divided into Blazhko and 
non-Blazhko RRLyrae stars, is [Fe/Hj^^^^^ = — 1.23dex. 
The difference in the mean metallicities between the two 
types (A[Fe/H] = 0.06 dex) is much smaller than the 
uncertainties associated with our method. We conclude 
that we would not detect a significant difference between 
the two types of RR Lyrae stars and therefore treat the 
stars independently of their type. Although star-to- star 
fluctuations in metallicity are observed on small scales, 
we find no evidence for a spatial metallicity gradient in 
the LMC. 



For the SMC we find a mean metallicity of 

[Fe/H]^^^^ = -1.42 dex on the [J95| scale based on 1831 
RRLyrae stars of typeab. This translate s to a v alue of 
[Fe/H]^^^^ = -1.70 dex on the scale of |ZW84ll and is 
in very good agreement with the spectroscopic deter- 
minations by Butler et al.^ {^^^ . '^) photometric 
metallicities by [Kapakos et al. T^Ol 1 ) . No metallicity 



gradient is observed in the old population of the SMC. 

In the LMC, several old globular clusters with ages 
greater than 10 Gyr are known. The metallicities of 



these clu sters span at least one dex in [Fe/H] (e.g, John- 



son et al.||2006| |Mucciarelh et al.||2010| ). This is smaller 



than the range of metallicities observed for Galactic 
halo globular clusters, but well within the range of what 
we found for the RRLyrae field stars, although most of 
the old LMC globulars with spectroscopic metallicity 
determinations are more metal-poor than the peak of 
our MDF of the RRLyrae stars. In the SMC only on e 
globular cluster, NGC 121, is known jGlatt et al.|2008b ) 
and h as a metallicity of —1.46 ± 0. 
scale (Da Costa & Hatzidimitriou 



_dex on the |ZW84 
1998D, weh within 



the range of metallicities of SMC held RRLyrae stars. 
The peak of our SMC RRLyrae MDF is slightly more 
metal-poor than the metallicity of NGC 121. The SMC 
in fact s hows a wide range of metallicities at a ny given 
age (e.g., |Kayser et al.||2007l [Glatt et al.||2008aD . 

Large metallicity spreads have also been found in 
most Local Group dwarf spheroidal (dSphs), usually 
exceeding one dex even when focusing on t heir old (red 



giant branch) populations (e.g., ^Grebel et al. 2003). 
For the old population of the LMC we find an intrinsic 
metallicity spread of a = 0.24 dex (in terms of one stan- 
dard deviation of the distribution), and a = 0.27 dex 
for the SMC. The full range of metallicities of the 
RRLyrae stars exceeds 1 dex in [Fe/H] and is difficult to 
constrain owing to the limits of the validity range of our 
calibration. As pointed out already, in neither of the 
two Clouds do we find a statistically significant trend in 
metallicity with position despite the spread. 

While we cannot assign specific ages to our individual 
RRLyrae stars, we suggest that this spread was caused 
by inhomogeneous, localized enrichment, which may 
hap pen on relatively short time scales as discussed 



by Marcolini et al. (2008). Detailed spectroscopic 
element abundance :atios measured for red giants in 
old LMC globular clusters indicate star formation from 
gas enriched by supernovae of Typ e II with prominent 
contributions from the r-process (jMucciarelli et aL 



2010). The a enhancement is comparable to Galactic 
halo stars (and thus inconsistent with the [a/Fe] ratios 
found in the old populations of many of the Galactic 
dSph and ultra- faint dSph stars; e.g., [Koch et aL] 
(|2008| and [Aden et al.] |2011). Jo hnson et al. (2006|) 
hnd differences between LMC globular cluster element 
abundance ratios both in comparison to dSphs and to 
Galactic globular clusters, emphasizing that the LMC 
experienced a different star formation history than the 
Galactic old populations. No such detailed, age-dateable 
information is available for the SMC thus far. 

If we qualitatively compare our RRLyrae MDFs to 
MDFs derived spe ctroscopically for G alactic dSph s (see, 
e.g.. Figure 13 in |Koch et al.] [20061 Fi gure 7 in [Koch 
et al.'2007a; Figure 12 in |Koch et al.|2QQ7 b, and Figure 1 
in ,Kirby et al., ,20TT] ) we find that the MDFs of the 
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Clouds do not show the steep drop-off at the metal-rich 
end found in some d Sphs. In the chemical evolution 



models for dSphs of Lanfranchi & Matteucci ([2004), 



such a steep decline is caused by gas outflows, which 
then do not seem to have played a defining role in the 
evolution of the old population of the Magellanic Clouds. 
This may not be surprising given the presumably much 
deeper potential wells of the Clouds. The shape of 
the MDFs of the MCs' old populations nevertheless 
resembles the shape of the MDFs of several G alactic 
dSphs that, according to |Kirby et al. 



well with the gas -infall models ot Lynden-Bel 



(20111 



and Pagel (19971), suggesting gas accretion from the 
intergalactic medium at early times. Note that since 
many of the Galactic dSphs contain predominantly 
old populations, we are comparing stars formed during 
similar epochs here. 

Although the MDFs of the RRLyrae stars in the 
MCs reveal a considerable metallicity spread, we do 
not see evidence for a large population of stars with 
[Fe/H] < — 2dex. We caution that the calibration of 
the Fourier decomposition method does not extend to 
low metallicities and that it has been suggested that 
this method overestimates metallicities at low [Fe/H] 



values (|Dekany fc Kovacs 2009 but see also Arellano 
1^'erro et al.||201l[ ). Additionally the number of KK Lyrae 
stars with low metallicities could be infiuenced by 
the shift of the instability strip to lower temperatures 
with lower metallicities. This could result in fewer 
horizontal branch s tars experiencing the RR Lyrae phase 
dYoon fc Leel|2002^ . However, very metal poor globular 
clusters (e.g. IN GC 2419 and NGC 7078 ) also contain a 
considerable fraction of RRLyrae stars ( [Catel an 2009). 
On the other hand, there are a number of metal-poor 
globular clusters both in the Milky Way and in the 
LMC that are dominated by blue horizontal branch 
morphologies a nd that essent ially do not contain any 
RRLyrae stars (Catelan 2009). Hence it is conceivable 
that the apparent lack of metal-poor RRLyrae stars is 
also caused by this evolutionary effect. Nonetheless, the 
apparent absence of a substantial, very metal-poor old 
population suggests that the early enrichment in the 
MCs proceeded rapidly and efficiently. The MCs show 
a similar lack of metal-deficient stars as the Galactic 
dSphs and the Galactic halo (also known as the "G 
dwarf problem"), suggesting that the majority of the 
old stars formed from pre-enriched gas. Spectroscopic 
studies of our RR Lyrae metal-poor candidates would be 
desirable to confirm their metallicities. 

The differences in the mean metallicities of the old 
populations in the LMC and SMC makes a common 
origin of these two galaxies from a single protogalactic 
fragment unlikely, unless such a putative progenitor 
would have be en separated into two ob jects very early 



(2010) proposed that 



on. Recently, Yang fc Hammer 
the LMC might have been ejected from M31 in a major 
merger event. The comparativel y high metallicity of 
old globular clusters in M31 (Caldwell et al.||2011 ) as 
compared to the properties of the MDF of the old LMC 
population seems to argue against this scenario. 
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